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1. INTRODUCTION 

It is generally believed that QCD with massless quarks undergoes a chiral 
phase transition (see jjj for a review). This leads to important observable 
signatures in the real world with two light quarks. The order parameter of the 
chiral phase transition is the chiral condensate (ipip), is directly related to the 
average spectral density of the Dirac operator Q. However, the eigenvalues 
of the Dirac operator fluctuate about their average position. The question we 
wish to address in this lecture is to what extent such fluctuations are universal. 
If that is the case, they do not depend on the full QCD dynamics and can be 
obtained from a much simpler chiral Random Matrix Theory (chRMT) with 
the global symmetries of the QCD partition function. 

This conjecture has its origin in the study of spectra of complex systems 
According to the Bohigas conjecture, spectral correlations of classically 
chaotic quantum systems are given by RMT. A first argument in favor of the 
universality in Dirac spectra came from the analysis of the finite volume QCD 
partition function. In particular, for box size L in the range 1/A Cl< 
(A is a typical hadronic scale and is the pion mass) we expect that the 
global symmetries determine its mass dependence j|. This implies that the 
fluctuations of Dirac spectra are constrained by Leutwyler-Smilga sum rules. 

Recently, it has become possible to obtain all eigenvalues of the lattice QCD 
Dirac operator on reasonably large lattices || . This makes a direct verification 



2 



J.J.M. Verbaarschot 



of the above conjecture possible. This is the main objective of this lecture. 

At nonzero chemical potential the QCD Dirac operator is nonhermitean with 
eigenvalues scattered in the complex plane. The possibility of a new type of 
universal behavior in this case will be discussed at the end of this lecture. 

2. THE DIRAC SPECTRUM 

The order parameter of the chiral phase transition, (ipip), is nonzero only below 
the critical temperature. As was shown in [Q (ipi/)) is directly related to the 
eigenvalue density of the QCD Dirac operator per unit four- volume 

|(^)| = ^. (1) 

It is elementary to derive this relation. The Euclidean Dirac operator for gauge 
field configuration is given by D = 7 M (<9 M + iA^). For Hermitean gamma 
matrices D is anti-hcrmitean with purely imaginary eigenvalues, D(f>k — iXkfik, 
and spectral density given by p(X) = J2k$(^ — ^k)- Because {75,!?} = 0, 
nonzero eigenvalues occur in pairs ±Afc. In terms of the eigenvalues of D the 
QCD partition function for Nf flavors of mass m can then be written as 

Z(m) = ([[(Xl + m 2 ) N f exp(-SV M )>, (2) 
k 

where the average (•) is over all gauge field configurations. 

The chiral condensate follows immediately from the partition function (Q), 

W)~5„JogZM = I(^^). (3) 

If we express the sum as an integral over the spectral density, and take the 
thermodynamic limit before the chiral limit so that we have many eigenvalues 
less than m we recover (|l|) (Notice the order of the limits.). 

An important consequence of the Bank-Casher formula (|l|) is that the eigen- 
values near zero virtuality are spaced as AA = l/p(0) = ir/Y,V. This should be 
contrasted with the eigenvalue spectrum of the non-interacting Dirac operator. 
Then p frcc (A) ~ ^A 3 which leads to an eigenvalue spacing of AA ~ l/V 1 ' 4 . 
Clearly, the presence of gauge fields leads to a strong modification of the spec- 
trum near zero virtuality. Strong interactions result in the coupling of many 
degrees of freedom leading to extended states and correlated eigenvalues. On 
the other hand, for uncorrelated eigenvalues, the eigenvalue distribution fac- 
torizes and we have p(A) ~ X 2Nf+1 , i.e. no breaking of chiral symmetry. 

Because the QCD Dirac spectrum is symmetric about zero, we have two 
different types of eigenvalue correlations: correlations in the bulk of the spec- 
trum and spectral correlations near zero virtuality. In the context of chiral 
symmetry we wish to study the spectral density near zero virtuality. Because 
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the eigenvalues are spaced as 1/£V it is natural to introduce the microscopic 



The dependence on the macroscopic variable £ has been eliminated and there- 
fore ps(u) is a perfect candidate for a universal function. 

3. SPECTRAL UNIVERSALITY 

Spectra for a wide range of complex quantum systems have been studied both 
experimentally and numerically (see Q for a review) . One basic observation has 
been that the scale of variations of the average spectral density and the scale 
of the spectral fluctuations separate. This allows us to unfold the spectrum, 
i.e. we rescale the spectrum in units of the local average level spacing. The 
fluctuations of the unfolded spectrum can be measured by suitable statistics. 
We will consider the nearest neighbor spacing distribution, P(S), the number 
variance, £2(71), and the As(n) statistic. The number variance is defined as 
the variance of the number of levels in a stretch of the spectrum that contains 
n levels on average, and As(n) is obtained by a smoothening of £2(71). 

These statistics can be obtained analytically for the invariant random ma- 
trix ensembles defined as ensembles of Hermitean matrices with independently 
distributed Gaussian matrix elements. Depending on the anti-unitary sym- 
metry, the matrix elements are real, complex or quaternion real. The corre- 
sponding Dyson index is given by (5 — 1, 2, and 4, respectively. The near- 
est neighbor spacing distribution is given by P(S) ~ S@ exp(— apS 2 ). The 
asymptotic behavior of £2(71) and A^n) is given by £2(71) ~ (2/tt 2 P) log(n) 
and A^n) ~ /3£2(n)/2. For uncorrelated eigenvalues one finds that P(S) — 
exp(— 5 1 ), £2(71) = n and A^n) — n/15. Characteristic features of random ma- 
trix correlations are level repulsion at short distances and a strong suppression 
of fluctuations at large distances. 

The main conclusion of numerous studies of eigenvalue spectra is that spec- 
tral correlations of a classically chaotic systems are given by RMT (Jt| || . 

4. CHIRAL RANDOM MATRIX THEORY 

In this section we will introduce an instanton liquid inspired RMT for the QCD 
partition function. In the spirit of the invariant random matrix ensembles we 
construct a model for the Dirac operator with the global symmetries of the QCD 
partition function as input, but otherwise Gaussian random matrix elements. 
The chRMT that obeys these conditions is defined by [|, [To| 



spectral density Q 




(4) 




/=i 
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and If is a n x m matrix with v = |n — m\ and iV = n + m. The matrix 
elements of W are either real (j3 = 1, chiral Gaussian Orthogonal Ensemble 
(chGOE)), complex ((3 — 2, chiral Gaussian Unitary Ensemble (chGUE)), or 
quaternion real ((3 = 4, chiral Gaussian Symplectic Ensemble (chGSE)). 

This model reproduces the following symmetries of the QCD partition func- 
tion: i) The Ua(1) symmetry. All nonzero eigenvalues of the random matrix 
Dirac operator occur in pairs ±A. ii) The topological structure of the QCD 
partition function. The Dirac matrix has exactly \v\ = \n — m\ zero eigenvalues. 
This identifies v as the topological sector of the model, in) The flavor symme- 
try is the same as in QCD. iv) The chiral symmetry is broken spontaneously 
with chiral condensate given by £ = limjv-*oc np(0) /N, (N is interpreted as the 
(dimensionless) volume of space time.) v) The anti-unitary symmetries. For 
fundamental fermions the matrix elements of the Dirac operator are complex 
for N c > 3 ((3 — 2) but can be chosen real for N c = 2 ((3 = 1). For adjoint 
fermions they can be arranged into real quaternions (f3 = 4). 

Note that spectral correlations of chRMT in the bulk of the spectrum are 
given by the invariant random matrix ensemble with the same value of (3. Both 
microscopic correlations near zero virtuality and in the bulk of the spectrum 
are stable against deformations of the ensemble. This has been shown by a 
variety of different arguments [jy] [l2| |]| |l5[ |l4| . 

Below we will discuss the microscopic spectral density. For N c = 3, Nt 
flavors and topological charge v it is given by || 

PS(U) = - (Ja( u ) - Ja+l(u)J a -l(u)) , (6) 

where a — Nf + v. The more complicated result for N c = 2 is given in |jl6|| . 

Together with the invariant random matrix ensembles, the chiral ensembles 
are part of a larger classification scheme. As pointed out in jnj , there is a one 
to one correspondence between random matrix theories and symmetric spaces. 



5. LATTICE QCD RESULTS 

Recently, Kalkreuter calculated all eigenvalues of the lattice Dirac operator 
both for Kogut-Susskind (KS) fermions and Wilson fermions for lattices as large 
as 12 4 . In the the case of SU(2) the anti-unitary symmetry of the KS and the 
Wilson Dirac operator is different jliL [HJ. For KS fermions the Dirac matrix 
can be arranged into real quaternions, whereas the Hermitean Wilson Dirac 
matrix 7 5 Z} Wllson can be chosen real. Therefore, we expect that the eigenvalue 



correlations are described by the GSE and the GOE, respectively 19 . In Fig. 
1 we show results for T, 2 {n), A 3 (n) and P(S). The results for KS fermions are 
for 4 dynamical flavors with ma — 0.05 on a I2 4 lattice. The results for Wilson 
fermion were obtained for two dynamical flavors on a 8 3 x 12 lattice. Other 
statistics are discussed in [G0|. 
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Recent lattice studies of the microscopic spectral density for quenched SU(2) 
show perfect agreement between the microscopic spectral density and random 
matrix theory These calculations were performed at moderately strong 

coupling. Result for couplings in the scaling regime are in progress. 

However, an alternative way to probe the Dirac spectrum is via the valence 
quark mass dependence of the condensate (22|, i.e. S(m) = J dXp(X)2m / (X 2 + 
m 2 ), for a fixed sea quark mass. In the mesoscopic range, S(m) can be obtained 
analytically from the microscopic spectral density (23| , 

^ = x(I a {x)K a (x) + I a+1 (x)K a ^x)), (7) 

where x — mVT, is the rescaled mass and a — Nf + v. In Fig. 2 we plot this 
ratio as a function of x for lattice data of two dynamical flavors with mass 
ma = 0.01 and N c = 3 on a 16 3 x 4 lattice. We observe that the lattice data for 
different values of (3 fall on a single curve. Moreover, in the mesoscopic range 
this curve coincides with the random matrix prediction for Nf = v = 0. 
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6. CHIRAL RANDOM MATRIX THEORY AT fi ^ 

At nonzero temperature T and chemical potential fi a schematic random matrix 
matrix model of the Dirac operator in (||) is given by pH p5[ " 



/ o iW + i n T + tl \ 
D - \ iw* + in T + m o J ' (8J 

where 57t = T®„ (2n+l)7iT. Below, we will discuss a model with fir absorbed 
in the random matrix and /i ^= 0. For the three values of (3 the eigenvalues of 
T> are scattered in the complex plane. 

In the quenched approximation the eigenvalue distribution for (3 — 2 was 
obtained analytically f27j from the Nf — ► limit of a partition function with the 
determinant replaced by its absolute value. The same analysis can be carried 
out for (3 = 1 and f3 = 4. In the normalization (5), it turns out that the average 
spectral density does not depend on (3 p6[ |. However, the fluctuations of the 
eigenvalues are (3 dependent. In particular, one finds a very different behavior 
close to the imaginary axis and not too large values of \i. This regime of almost 



Hermitean random matrices was first identified by Fyodorov et al. |28 . For 
(3 = 1 we observe an accumulation of purely imaginary eigenvalues, whereas for 
(3 = 4 we find a depletion of eigenvalues in this domain. These results explain 
quenched instanton liquid calculations Ej| and lattice QCD simulations with 
Kogut-Susskind fermions (30| at [i ^ (both for SU(2)), respectively. 



7. CONCLUSIONS 



We have shown that microscopic correlations of the QCD Dirac spectrum can 
be explained by RMT and have obtained an analytical understanding of the 
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distribution of the eigenvalues near zero. An extension of this model to nonzero 
chemical potential explains previously obtained lattice QCD results. 
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